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\ Abstract. We present an extension of the theory known as Lin's method to 

i-^h ■ heteroclinic chains that connect hyperbolic equilibria and hyperbolic periodic orbits. 

d . Based on the construction of a so-called Lin orbit, that is, a sequence of continuous 

partial orbits that only have jumps in a certain prescribed linear subspace, estimates 
for these jumps are derived. We use the jump estimates to discuss bifurcation equations 
for homoclinic orbits near heteroclinic cycles between an equilibrium and a periodic 
> ' orbit (EtoP cycles). 

(N 
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; AMS classification scheme numbers: 37C29, 37G25, 34C23, 34C60 
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Q\ , 1. Introduction 

O 

Connecting cycles involving hyperbolic equilibria and hyperbolic periodic orbits play an 
^ ■ important role in many applications, cf. [12] and references therein. The bifurcation 
■ analysis in the vicinity of such a connecting cycle is crucial for the understanding of 
the system's behavior. In this respect, both the theoretical bifurcation analysis and 
numerical implementations are of high interest in current research. 

Lin's method has proved to be an appropriate tool for discovering recurrent 
dynamics near a given cycle. The method dates from [II] . where heteroclinic chains 
consisting of hyperbolic fixed points, all having the same index (dimension of the 
unstable manifold), and heteroclinic orbits connecting them are considered. The basic 
idea of Lin's method is to construct discontinuous orbits with well defined discontinuities 
(jumps), so-called Lin orbits, near the original cycle. By 'making these jumps zero' one 
finally finds real orbits staying for all time close to the cycle under consideration. In 
1993 Sandstede [19] gained jump estimates which allow an effective discussion of the 
bifurcation equations. For a survey of the many applications and several extensions of 
the method we refer to H5l. 
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In this paper we present an extension of the theory of Lin's method to arbitrary 
heteroclinic chains connecting hyperbolic equilibria and hyperbolic periodic orbits. 
Related problems have been studied in [5j [16j [TTj and [IE]; for the numerical 
implementation of these ideas we refer to |12j . 

In the presence of periodic orbits the construction of Lin orbits is much more 
involved, because the dynamics near the periodic orbit has to be incorporated. The 
handling of the flow near the periodic is the main difference between the approaches 
in jT6l [IT] and p2]. Based on the ideas in [18], we construct certain partial 
(discontinuous) orbits running between Poincare sections of two consecutive periodic 
orbits of the given chain. Then the dynamics near each periodic orbit is described by 
means of the corresponding Poincare map. Finally, the different orbits are coupled in 
the Poincare section in each case. This approach allows to apply immediately results 
from Lin's method for discrete systems |10j . 

We consider a family of ODE 



For a particular parameter value, say A = 0, we assume that the system has a heteroclinic 
chain consisting of hyperbolic periodic orbits 7; and heteroclinic orbits qi connecting 7, 
and 7i+i. Here we explicitely admit that the minimal period of either of these periodic 
orbits may be zero, meaning that either of these orbits may be an equilibrium. We want 
to note that, for instance, a heteroclinic cycle between an equilibrium and a periodic 
orbit can be considered as such a heteroclinic chain. In this case the chain consists of 
copies of the cycle under consideration which are stringed together. 

We refer to a segment 7$ U qi U 7^+1 of the given chain as a short heteroclinic chain 
segment. Near q^ we construct a discontinuous orbit X-i satisfying certain boundary 
conditions (B~) near and (B^ +1 ) near 7i+i. There the discontinuity is a well defined 
jump Sj near ^ (0) . Those orbits we call short Lin orbit segments. 

It can be shown that arbitrarily many consecutive short Lin orbit segments can be 
linked together to a Lin orbit close to the original chain, see [18] for chains related to 
heteroclinic cycles connecting one equilibrium and one periodic orbit. In the present 
paper we confine ourselves to linking two consecutive short Lin orbit segments X\ and 
X r related to 7; U qi U 7 and 7 U q r U 7 r to a long Lin orbit segment with boundary 
conditions (Bf) and (-B+). Apart from the fact that this procedure reveals the basic 
idea for linking arbitrarily many consecutive short Lin orbit segments, it is eligible for 
consideration in its own right. So it suffices to consider long Lin orbit segments for the 
detection of 1-homoclinic orbits near a heteroclinic cycle connecting two periodic orbits. 
Here, 1-homoclinic orbits are characterized by only one large excursion before returning 
to their starting point. 

If 7 is an equilibrium, the existence proof of long Lin orbit segments runs to large 
extent parallel to 'classical constructions' of Lin's method, see [18]. For that reason we 
consider only the case that 7 is a periodic orbit with nonzero minimal period. Roughly 
speaking, the orbits Xi and X r are linked via an orbit x that defines the behavior of 
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the newly generated orbit along 7. We construct x = x(y) as a suspension of a certain 
orbit y of an appropriate Poincare map. In this process the boundary conditions (Bf) 
and (B+) remain untouched. 

Finally we give estimates of the jumps and S n which allow us to discuss the 
bifurcation equations Sj = and H r = for detecting actual orbit segments near the 
given long orbit segment 7; U qi U 7 U q r U 7,.. 

We apply our results to study homoclinic orbits near a heteroclinic cycle connecting 
a hyperbolic equilibrium E and a hyperbolic periodic orbit P (with nonzero minimal 
period), an EtoP cycle for short. Here we only consider 1-homoclinic orbits to the 
equilibrium. Those orbits may differ considerably in their length of stay near P. This 
length correlates to the number v of rotations the homoclinic orbit performs along P 
or, in the above notation, it correlates to the length v of the orbit y. 

Indeed, in numerical computations it has been observed that the homoclinic orbits 
for different v all lie on the same continuation curve. Moreover, this continuation curve 
shows a certain snaking behavior and accumulates on a curve segment related to the 
existence of the primary EtoP cycle, cf. figure [5] panel (a). The addressed snaking 
behavior of a system with reinjection was revealed numerically in [Tl] and [12] . 

In particular, we explain two local phenomena appearing in this global snaking 
scenario. First we consider a co dimension-one EtoP cycle. Apart from E and P this 
cycle consists of a robust heteroclinic orbit q\ connecting E to P, and a codimension- 
one heteroclinic orbit q r connecting P to E. Further, the dimensions of the unstable 
manifold of P and the stable manifold of E add up to the space dimension. 

Let A be the one-dimensional parameter unfolding q r and, hence, unfolding the 
entire cycle. In that unfolding we describe the accumulation of homoclinic orbits at the 
primary heteroclinic cycle. More precisely, we prove that there is a sequence (\ v ) u m 
tending to zero such that for each \ v there is a homoclinic orbit to the equilibrium, while 
for A = the heteroclinic cycle exists. Moreover, with increasing v the corresponding 
homoclinic orbits stay longer near the periodic orbit, performing an increasing number 
of rotations along the periodic orbit. 

In a second scenario we assume that W U (E) and W S (P) do no longer intersect 
transversally but have a quadratic tangency — still we assume that q r is of codimension 
one as described above. Then the entire EtoP cycle is of codimension two. Let 
the parameter Xu r unfold the orbits qu r , and assume that the EtoP cycle exists for 
A = (A;, A r ) = 0. Then, in the neighborhood of A = 0, we find a sequence of curves 
K y in the A-plane for which a homoclinic orbit to the equilibrium exists. As above 
with increasing v the corresponding homoclinic orbits stay longer and longer near the 
periodic orbit, performing an increasing number of rotations along the periodic orbit. 
For each v the curve k u has a turning point X u tending to zero as v tends to infinity. 
This explains the curve progression of hi in a small neighborhood of c\ n to? cf- figure [5] 

The paper is organized as follows. In section [2] we develop Lin's method for short 
heteroclinic chain segments. The main theorem in this respect is theorem 12.21 which 
states the existence of short Lin orbit segments. Corollary 12.81 extends theorem 12.21 to 
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boundary conditions that enforce that the short Lin orbit segment 'starts' in the unstable 
manifold of 7$. In lemma [2. 101 we give an estimate of the jump function, which we extend 
in corollary 12.111 to the situation of corollary 12.81 In section [3] we describe the coupling 
of two consecutive short Lin orbit segments Xi and X r to a long Lin orbit segment. 
Their existence is stated in theorem 13.11 and corollary 13.21 extends this assertion on Xi 
and X r asymptotic to 7; and j r , respectively. The corresponding estimates of the jumps 
are given in lemma 13. Ill and corollary 13. 121 In section H] we consider 1-homoclinic orbits 
near EtoP cycles. The corollaries 14.31 and 14.41 explain the accumulation of 1-homoclinic 
orbits, and corollary 14.71 describes the accumulation of vertices of continuation curves 
that are obtained by unfolding a ray (in parameter space) that is covered twice. 

2. Lin's method for short heteroclinic chains 

Consider the ODE (00). Throughout this section we assume that for A = there is a 
short heteroclinic chain segment 7~ U q U 7 + with hyperbolic periodic orbits 7" and 
7 + . We explicitly admit that the minimal period T~/ + of either of them may be 
zero, meaning that 7" and/or 7 + may be hyperbolic equilibria. Let (7*) denote 
the stable/unstable manifolds of 7 ± , and we use the short notation iy s/ ' M (7 ± ) for the 
corresponding manifolds at A = 0. Further, T q W s / u denotes the tangent space of the 
corresponding manifold at q. 

We introduce subspaces W + , W~ and U as follows: 

(T m w u (r) n T m w s ( 7 + )) ®w- = T m w u (r), 

(T q{0) W u (j-) n T g(0) iy s ( 7 + )) © W + = T 9(0) ^( 7 + ) and 
span{/(g(0), 0)} © U = T^W^) R T, (0) W s ( 7 +). 

In other words, the linear spaces W~ and W + are contained in the tangent spaces of the 
unstable and stable manifolds of 7" and 7 + at q(0), but do not contain their common 
directions, which are collected in span{/(g(0), 0)} © U. 
Using a scalar product (■, ■) in IR n we define 

Z := {W + © W~ © U © span{/(g(0), 0)}) x (2) 

and 

Y = W + © W~ © U © Z, (3) 

and we denote the projection onto U in accordance with the direct sum decomposition 
([3]) by P u . Note that either of the involved spaces W ± , U and Z may be trivial. Finally, 
we define a cross-section E of q as 

E := q(0) + Y. 

Our goal is to construct 'discontinuous orbits' near q that satisfy certain boundary 
conditions (B~) and {B + ) near 7" and 7 + . Actually those orbits consist of two orbit 
segments where the end point of the first and the starting point of the second one are 
in E and their difference is in Z, which is reflected in the boundary condition (J). 
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In a first step we prove the existence of orbit segments that lie in the unstable and 
stable manifolds of 7" and 7 + , respectively, and that satisfy certain jump conditions 
in E. 

Theorem 2.1. There is a constant c > such that for |A| < c and u G U , \u\ < c, 

there is a unique pair of solutions (q~(u, A), q + (u, A)) of (CP that satisfy 

(i) q + (u,\)(0) e W^+), q-(u,\)(0) E M^(t"), 
(i%) g + ( M ,A)(0),g-( M ,A)(0)GS ; 

(in) P u (q + (u, A)(0) - q(0)) = P u (q~(u, A)(0) - q(0)) = u and 
(iv) q+(u,\)(0)-q-(u,\)(0)eZ. 

Figures [1] and [2] give a graphical interpretation of theorem 12 .11 the proof is given in 
section 12. 11 




Figure 1. Sketch of the situation described in theorem 12.11 showing the orbits 
q~ C W£(j~) and q + C W^(j + ). Within the cross-section E, the two orbits have 
a jump in the direction Z. Note that 7" and 7 + arc depicted as periodic orbits, but 
cither of them may be an equilibrium. The dotted connection q is present for A = 0. 




Figure 2. The situation inside E, depicted are the traces of q~/ + {u, A), W^-f^), 
W^(j + ) and the direction U consisting of the common tangent directions (restricted 
to E). The depicted situation corresponds to a 'quadratic tangency' of W u (^~) and 
W s {j+). 
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In the next step, we perturb the solutions given by theorem 12.11 to construct 
solutions that stay near the connecting orbit q and satisfy projection boundary 
conditions near 7" and 7 + . Moreover, these solutions are also allowed to have a jump 
in the direction Z (within E). 

To formulate the boundary conditions we define projections P ± (u, X)(t) by 



imP + (w,A)(0) =T C 



q+(u,X)(0) 



w s x (r 



kerP + (u,A)(0) = W~ © Z and 

P + (u,X)(t) := $+(t,0)(P + (u,A)(0))$ + (0,t), t G 
and analogously 



(4) 



imP-(u,A)(0) =T C 



9-(u,A)(0) 



(5) 



kerP~(n, A)(0) = W + © Z and 

p-(u,X)(t) := $-(t,o)(p-(u,A)(0))$-(0,t), t g R". 

Here $ ± (-,-) = $ ± (m, A)(-, •) denotes the transition matrix of the variational equation 
along g ± (u, A)(-). 

Throughout we denote by |-| the absolute value of a number or the the Euclidian 
norm of an n-tuple. For elements a = (a~,a + ) G M n x M. n we define ||a|| := 
max{|a~|, |a + |}. 

Theorem 2.2. Fix uj + > 0. There is a constant c > such that for |A| < c, u G U , 

\u\ < c, and a = (a~,a + ) G M n x W 1 , \\a\\ < c, there is a unique pair of solutions 
(x~,x + ) of (TJP that satisfies 

(J) x~/ + (a, u, A)(0) G E, x~(a, w, A) (0) — x + (a, u, A) (0) G 
(£-) (id-P-(«,A)(-a;-)) (ar _ (a,u, A)(-w~) - g~(u,A)(-a;-) - a - ) = 0, 
(P+) (^-P + (M,A)(cu + ))(x + (a,M,A)(cj + ) -q+(u,X){u} + ) - a+) =0. 
Moreover, (x~,x + ) = (a, w, A) depends smoothly on (a,u,X). 




im (id- P-(u,X)(-u-)) 
\ 

\ 



im (id-P+(w,A)(a;+)) 
/ 

/ 



Figure 3. Sketch of a short Lin orbit segment (x , x + ) near a short heteroclinic chain 
segment 7~ U g U 7+. 



In other words, theorem 12.21 states the existence of short Lin orbits segments 
X = with boundary conditions (B~) and (B + ) for sufficiently small (a,u, A), 

for a sketch of this situation see figure [3j Note that X also depends on u~,uj + > 0, 
which we suppress from our notation in this section. The proof is given in section 12.21 
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2.1. Orbits in the stable /unstable manifolds — the proof of theorem \2.1\ 

From a geometrical point of view, the statement of theorem 12.11 is rather clear. Hence 
we give a proof which exploits the geometry of the traces of the involved manifolds in S, 
see figure [2j By Bx{x,r) we denote a closed ball in X centered at x with radius r. If 
the space X is clear from the context we will also write B(x,r) for short. 

Proof of theorem \2.1\ Using the direct sum decomposition ([3]) we find the following 
representations of the traces in £ of the stable/unstable manifolds of 7"/+ locally around 
q(0): For e sufficiently small and for w s G W^a(7 + ) ^ £ H B(q(0),e) there are smooth 
functions w~ : W + x U x R m -> W~ with tiT (0,0,0) = 0, L>i?ir(0, 0, 0) = and 
z + : W + x U x R m -> Z such that 

w s = q(0) +w + + w~(w + , w + , A) + z + (w + , u + , A) + u + . 

Similarly for w u G W^( , ~f~) PI £ H B(g(0), e) there are smooth functions w + : x U x 
_> w ith ^+(0,0,0) = 0, DitZ)+(0, 0, 0) = and z" : x(7x M m -> Z such 

that 

w" = g(0) + u - , A) + u; - + u~ , A) + 

The demand that q + {u, A)(0) — q~(u, A)(0) G Z results in u~ = u + =: u and 
w + = w + (w~, u, A), 
w~ = w~(w + , u, A), 

which then can be solved for (w + ,w~) = (w + (u, A), w~(u, A)) around (u, A) = (0,0). 
Now we get A) as the solution of the initial the value problem 

x = f(x, A) 

x(0) = g(0) + w + (u, A) + A), w, A) + A), w, A) + m, 

and q + (u, A) is the solution of the initial the value problem 

% =f(x,X) 

x(0) = q(0) + w + (w~(u, A), u, A) + w~(u, A) + z~(w~(u, X),u, A) + u. 

□ 

S.^. Short Lin orbit segments — the proof of theorem \2.2\ 

In this section we give a detailed proof of theorem 12.21 This proof is based on the 
ideas of Lin's method, but used in a slightly different way. The main difference to the 
'classical' proof of Lin's method is that we keep the boundary conditions near 7" and 
7 + as linear projection conditions, while finding solutions of the full nonlinear system 
that additionally satisfy certain jump conditions. 

We start with a lemma that provides some properties of the projections P ± (m, A)(-), 
as introduced in (Tj0) and ([5]), that are used in the proofs throughout this section. 
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Lemma 2.3. There are projections P+(tt, A)(-) and P c + (w, A)(-) such that 

P+(u, X)(t) = P+(u, X)(t) + P+(u, A)(i) for t G IR + . 
The projections P^ c (u, A)(-) satisfy the following: 

(i) 9(t, r)P+ c (r) = P+Jt)Q{t, r) Vt, r G M+ ; 

(ii) there are constants K > 0, 5 s , 5 U > 5 C > such that 

\${t,T)P+{r)\ < Ke- ss{t - T) ,t>r, 
\<S>(t,T)P+(r)\ < Ke &c{t - T) ,t>r, 
|$(t,r)P+(r)| < Ke~ sc{t - r) ,r>t, 
\<S>(t,T)(id-P + (r))\ < Ke- &u{T - l \T> t. 

This lemma follows immediately from the fact that the variational equation along 
the solutions q + (u, A) has an exponential trichotomy on M + [3j [9], E]. Note that the 
exponents S s / C / U are determined by the eigenvalues/Floquet exponents of 7 + . Since 7 + 
is a hyperbolic periodic orbit, we have 5 C = [3]. We want to note explicitely that 
the images of P+(u, \)(t) are well-determined — these are the tangent spaces of the 
strong stable fiber of 7 + at q + (u, A)(£). Also note that if 7 + is a hyperbolic equilibrium, 
the variational equation along q + (u, A) has in fact an exponential dichotomy [7J, i.e. 
P+(u, A) = 0, and imP s +( M , X)(t) = T q+(u>m W^+). 

For the projection P~(u,X) an analogous lemma holds, exploiting that the 
variational equation along the solutions q~(u, A) has an exponential trichotomy on R~: 

P-(«,A)(*) =P~(n,A)(t)+P~(u,A)(t) forteR" 

where P~(u, X)(t) projects on the tangent space of the strong unstable fiber at 
q-(u,X)(t). 

To prove theorem 12.21 we consider small perturbations of the solutions g~/ + (w, A): 

x-(t):=q-(u,\){t) + v-(t), teR" 

x + (t) :=q + (u,X)(t) + v + (t), tER + . 

The perturbation terms v~ and v + are solutions of 

v~ = D 1 f(q-(u,X)(t),X)v~ + h~(t,v~,u, A), 

i>+ = DJiq+iu, A)(t), A)v+ + /i+(t, u+, u, A), 

where h±(t, v, u, A) := /(g±(«, A)(t) + v, A) -/(</*(«, A)(t), X)-D 1 f(q±(u, X)(t), A). The 
boundary conditions (J) and (P ± ) for x 1 * 1 yield boundary conditions for v^: 

(J v ) v ± (0) G W~ © W + © Z, u-(0) - u+(0) G Z, 

(id - P"(u, A)(-w")) ^(-W) = (id - P-(u, X)(-u-)) a~, 
(P+) (id-P+(n,A)(^+))^+(w + ) = (id-P + (u,A)(cu+))a+. 

In a first approximation of (J7J), we replace the functions /i~ and /i + by functions g~ 
and g + only depending on t: 

v- = Dtf(q-(u, X)(t), X)v~ + g-(t), 

v + = D 1 f(q + (u,X)(t)A)v + + g + (t). 
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For given uj~ and u + we write oj := (uj~,lj + ), and we introduce the space V w of 
pairs of continuous functions as 

V u := { (v~, v + ) : v- G C ( [-0J-, 0] , R") and v + e C { [0, u + ] , R n ) } . 

We equip V u with the norm || (v~ ,v + ) \\ := max{||u~||, ||w + ||}, where denotes the 

supremum norm. 

We actually perform the proof of theorem 12.21 in two steps. First we consider the 
'linearized' equation (jSJ) (linearized in the sense that g does not depend on v^) and 
show that there are unique solutions satisfying boundary conditions (J v ), (B~) and 
(By); see Lemma [2.41 below. Of course v^ 1 depend (among others) on g^. In the next 
step we replace the function g ± in these dependencies by h ± . This gives a fixed point 
equation (see (TTB|) below) that is equivalent to (JTj) with boundary conditions (J v ), {B~) 
and (J3+). 

Lemma 2.4. Let u and X be in accordance with theorem \2.1[ and let uj be fixed. Then, 
for given a = (a~,a + ) G M. n x R n and g = (g~,g + ) G V w there is a unique pair of 
solutions v = (v~,i) + ) G V w , v = v(g,a,u, X), of (Uj) that satisfy boundary conditions 
(J v ), (B~) and (B+). 

For fixed u and X the pair of solutions v depends linearly on (g,a), and it depends 
smoothly on (g,a,u,X). Moreover, there are constants C a , C g > such that 

\\v(g,a,u,X)\\ < C a \\a\\ + C g \\g\\. (9) 

The constant C a is uniform in uj + and uj~ , while C g is uniform in lj^ only if 7 ± is an 
equilibrium. 

Proof. For a shorter notation in this proof we omit the dependencies of <3> and P^ on u 
and A. The variation of constants formula gives 

ir(*) = $-(i,0)iT(0)+ / $-(t,T)g-(r)dT, 

J \ (10) 

v +(t) = $ + (t,0)t; + (0) + / <S> + {t,r)g + (T)dT 

Jo 

as solutions of (jHJ), which can be transformed to 

(id - P~(0)) iT(0) = $"(0, -uj') (id - p-(-ur)) v~(-u~) 

+ / <^-(0,r) (id- p-(r))g~(r)dr, 

J — U) 

(id-P + (0))i; + (0) = $+(0,^+) (id-P + (u + ))v + (u + ) 



f $ + (0, r) (id -P + (r))g + (r)dr. 
Jo 



Heteroclinic chains involving periodic orbits 
Thus 

(id - P~(0)) u _ (0) = $"(0, -a; - ) (id - p-(-uT)) cT 

+ T $"(0, r) (id -P-(r))^-(r)dr, 

J — U)~ 

(id-P + (0))*;+(0) = $+(0,cu+) (id-P+(cu + ))a + 

-/ $ + (0, r) (id- P + (r)) £ + (r)dr. 

We decompose v~(0),v + (0) in accordance with the boundary condition {J v 
v~(0) = w~ + w + + z~ , 
v+(0) = w ~ +w + + z + , 

where w~ G W~ , w + G W + and z^ 1 G Z. 

With that the left-hand side of ( fTTj) can be considered as a linear mapping 

L : lf + x If" x Z x 2 -> (W+ © Z) x (If" © Z) 
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11^ 



(w + + z~) 
(w~ + z + ) 



which is invertible. Hence we can solve (fTTj) for v + (0), v~(0) linearly depending on (g, a). 
Incorporating the dependence on (it, A) finally gives a solution v = v(g,a,u, A) of (|HJ) 
satisfying the boundary conditions (Jj,), (-B~) and (Pj~). Note that v depends linearly 
on (g, a), and smoothly on (g, a, u, A). 

To prove estimate (Q, we decompose v + by means of the projection P + : 

+ (t) = (id - P + (t))v + (t) + P + (t){) + (t). 



Thus we have 



\v + {t)\ < \ {id - P + {t))v + {t) \ + \P + (t)v + {t)\ . 



(12) 
to derive 



We use the variation of constants formula and the estimates of lemma 
an estimate for the second term of (1121) : 

\p+(t)v+(t)\ = P + (t) (V(t,0){) + (0) + j\ + {t,T)g + {T)dT 

< K(e- 5H + e SH ) \v + (0)\+M\\g + \\ 

< K{e- 6H + 1) \v + (0)\+ M\\g + \\. 

The constants 5 s , 5 C and K are the corresponding constants of the exponential 
trichotomy (5 s > 5 C = 0; see lemma 12.31 and the remarks above). Note that if 7 + 
is not an equilibrium, the constant M depends on u + , in fact M — > oo as u + — ► oo in 
the same order as w + . 

We estimate |i) + (0)| by applying L^ 1 to (fTTj) and using lemma [2731 once again: 



|£ + (0)| < \\L^\\K (|a + | + \a~\) + M\\(g + ,g- 



(13) 



Heteroclinic chains involving periodic orbits 1 1 

Here, the constant M is uniform in u + and oj~ . 
Thus we have 

\P + (t)v + (t)\ <Ci, a ||a|| + C ljfl ||0||). 

The constant C\ %a is uniform in u + and to~ , while Ci i9 (u}) tends to infinity in the same 
order as ||u>||. 

For the first term of the right hand side of (I12jl we use 
(id - P+(t))v + (t) w + )(id - P + (uj + ))a + 

-£ $ + (t,T)(id-P + (T))g + (T)dT 

and thus finally get 

\(id- P + (t))v + (t)\ <C 2 {\a + \ + \\g\\). 

Note that Ci is uniform in u + (and does not depend on ou~). 
Summarizing, there are constants C£ and C£ such that 

\\v + (g,a,u,X)\\<C+\\a\\+C+\\g\\, 

where is uniform in u and C£(u}) tends to infinity in the same order as 
Proceeding with v~ in a similar way yields estimate (Q. □ 

Lemma 2.5. Let the assumption of lemma \2.J\ hold. We define functions 
a>i(9, a, u, A) := P+(u, X)(uj + )v + (g, a, u, X)(uj + ), 
al(g,a,u,X) := P~ (u, X)(-u~)v~ (g, a,u, X)(-u~). 

There are constants 5 s , 5 U > and C > such that 

I at I < C(e- ssu>+ \\a\\ + \\g\\) , la^l < C (V* 5 "^ \\a\\ + \\g\\) . (14) 



For the derivatives of a~[_ and a ± the following estimates hold: 

\\D 2 a+(g, a, u, X) \\ < Ce^ + , \\D 2 al(g, a, u, X) || < Ce' 6 ^ (15) 

and 

\\Diat(g,a,u,X)\\, WD^fa, a, u, A)|| < (7. (16) 
Proof. For estimate 014|) we use lemma 12.31 again: 

|a+(^a,M,A)| = P+(cu+)(^ + (cu + ,0)v+(0)+ £ ^ + {u + ,T)g + {r)dr) 

< Ke- 5Su+ \v + (0) \ +M\\g + \\ 

<C(e- ss " + \\a\\ + \\g\\). 

For the derivative we note that the dependencies of v ± on (g, a) are linear. 
Hence, there are linear operators depending on u and A such that v ± (g,a,u, X) = 
L ± {u,X){g,a) = Z ± (m, A)((?, 0) + Z ± (m, A)(0, a). Due to their definition, also depend 
linearly on (g,a). Thus the estimates (fT5l) and (1T61) follow. □ 
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Now we look for solutions of (0) satisfying the boundary conditions («/„), (B~) and 
{B£). For that purpose, in v(g,a,u, A) we formally replace the function g by /i. To 
perform this substitution we define the Nemytskii operator H := (H~,H + ): 

H : V u x U x R m -> , . 

(v,«,A) ~ (tf-(^u,A),tf+( W +,u,A)), 1 j 

where 

£T (v, u, A)(-) := h~ (-, «(•), «, A) , u, A)(-) := fc+ (-, u(-), u, A) . 

In [19] it is verified that H has the stated mapping properties and that H is smooth. 



Summarizing, we find that a function v solves the boundary value problem (([7]), 
(Jo), (B~), (-B+)) if and only if it satisfies the following fixed point equation in V u : 

v = v(H(v, u, A), a, u, A) =: F(t>, a, u, A). (18) 

Note that 

F:V u x (M n xr)x[/xl m ^ VL,. 

Lemma 2.6. Fis some a;. There are functions e,c,c,Q : R — > M + sitc/i £/ia£ for 
all K > 1 £/ie /ixed pomi problem fT8\) has a unique solution v = (v~,v + ) G 14,, 
w = v(a,u,X), in an e(K) -neighborhood of G V^, provided that |A| , |it| < c(if) ; and 
||a|| < c(K). The solution v depends smoothly on (a,u, A). 

Proof. We use the Banach fixed point theorem to prove existence and uniqueness of v. 
First we show that there is an F-invariant closed ball B(0,e) C V w and then that F is 
a contraction on £>(0,e) with respect to v. 

Let C a , C g and C be the constants in accordance with lemma 12.41 and lemma 12. 5[ 
respectively. Define C := max{l, C, C a , C g }. Then, according to (jHJ), 

||F(u,a,u,A)|| < C(\\a\\ + \\H(v,u,\)\\). (19) 

We start with an estimate for \\H(v , it, A)||. From the definition of h we see that 
H(0, 0, 0) = and thus we can use the mean value theorem to get the estimate 

\\H(v,u,X)\\ < [ ||£>i#(s(v,u,A))||ds||v|| 

\\D2H(s(v, u, X))\\ds \u\ 



+ I \\D 3 H(s(v,u,X))\\ds\X\. (20) 
Jo 



Applying the mean value theorem to D\H we find that there is an appropriate constant 
D > such that with 

m := ^ (21) 

the following holds: If |A| , \\v\\, \u\ < e(K) then, since DiH(0, 0, 0) = 0, 

||A#(^,A)||<_^. (22) 
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Further, there is a constant E such that for all v, u and A taken from some neighborhood 
of the origin 

/ \\D 2 H(s(v,u,X))\\ds< E, [ \\D 3 H(s(v,u,X))\\ds < E. 
Jo Jo 

By means of e we further define 

c(K) - ' {K) c - (23) 

Therefore, we find for \\v\\ < e, |A| , |it| < c and K > 1 (recall C > 1) 

\\H{v,u,\)\\<——< 



7K 2 C 2 ~ 7 KG 

Finally, estimate (fT9l) yields that for all |A| , |u| < c and ||a|| < c the mapping 
F(-, a, u, A) leaves the closed ball B(0, e) C V w invariant. 

Moreover, due to the linear dependence of v on (H, a) and the estimates ([9]) and 
( 122|) . we have 

||D 1 F(T;,a,«,A)|| < WD^H, a,u, A)|| ■ ||£>i#K a,«, A)|| < < ^. 

Thus, F is a contraction on -8(0, e) and the existence and uniqueness part of the lemma 
follows immediately from the Banach fixed point theorem. 

Applying the implicit function theorem at a solution point of ( TT8l) provides the 
smooth dependence of v on (a, u, A). □ 

Lemma 2.7. Choose K > 1 and (a,u,X) in accordance with lemma \27b\ We define 
functions 

a^(a, u, A) := P*(u, X)(u + )v + (a, u, A)(cj + ), 

aj(a, u, A) := F~(w, A)(— u;~)u~(a, it, A)(— c<j~). 

There are constants C,5 S ,5 U > 0, that do not depend on K , such that with c according 
to (G2P 

\a+(a,u,X)\ <C||a||e^ + + ^, 

^ (24) 

\al(a,u,X)\ < C\\a\\e- 5U "~ + 

and 

\D ia +(a,u,X) \ < Ce"" 5 ^ + 

\ (25) 

\D ia l{a,u,X)\ < Ce~ 5 ^ + — . 

K 

Proof. From (|14p and the estimates in the proof of lemma 12.61 we get 
\al(a,u,X)\ < C (||a||e^ + + \\H(v, u, X)\\) 
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Next we estimate Dia~\_(a, u, A) = ^a\_(H{v(a, u, X),u, A), a, u, A). Hence 

\Dia\(a,u, A) | < \Dia]_{H, a, u, A) | • ||-Diif(t> , w, A)|| ■ ||Di?;(a, u, A)|| 

+ |D 2 ai(-H", a, w, A) | . 
Because v (a, u, A) = v(H(v(a, u, A), u, A), a, u, A) we get 

Z>lU(a,U,A) < rrrr. 

1 — ||-Dit> (ii, a, m, A) || • \\DiH{v, u, X) \\ 

Recall that \\D 2 v(H,a,u,X)\\ < C, \\Div(H,a,u,X)\\ < C and \\D 1 H(v, u, A) || < ^jtc, 
hence 

ii n ( \\\\* 7KC 
\\Div(a, u, X)\\ < 



7K-1 

Together with (fT5l) and (|T6l) we finally get 

|-DiOx(o,«,A)| < + < i + Ce"- + . 

With similar computations for a J, the estimates of the lemma follow. □ 

Recapitulating, we want to note that we find solutions according to theorem 12.21 by 
inserting the solutions v(a, u, A) into the representation (jBJ). 

The statement of theorem 12.21 remains true for = oo' in the following sense: 

Corollary 2.8. Fix u + . There is a constant c > such that for |A| < c, u £ U, \u\ < c, 
and a + £ W 1 , \a + \ < c, there is a unique pair of solutions (x~, x + ) of (TJp that satisfy 

(J) x"/+(a + ,u,A)(0) £ £, x-(a+,u,\)(0) £ W^~), 
x~(a + , u, A)(0) — x + (a + , u, A)(0) £ Z and 

(B+) {id - P + (u, X)(u+)) (z+(a+ u, X)(u+) - q+(u, X)(u+) - a+) = 0. 

Proof. Basically the statement follows by setting (id — P~(u, A)(— uj~)) a~ = in 
theorem 12. 2t Let uj~ be any value in accordance with theorem I2.2L Then, due 
to (id — P~(u, A)(— uj~)) a~ = 0, it follows that v~(— uj~) £ imP~(u, A)(— u~) = 
Tq-( u ,\)(-u>-)Wx( , y~), cf. Lemma 12.41 Assuming that W^{^~) is flat around 
q~(u, A)(— u)~), meaning that locally around q~(u, A)(— uo~) the unstable manifold 
W"(7~) and A)(-a;~) + T ? - (UiA) (_<,,- (7") coincide, we find q~(u,X)(-u~) + 

fi-(-a,-) e W^iT)- 

Solving fixed point equation (JT5J) with that particular 0, we find that 
x~(m, A)(— uj~) = q~(u, X)(—uj~) + v~(— u~) £ W% (7~); compare also (E]). Hence x~ lies 
in the unstable manifold of 7~. □ 

Remark 2.9. In the same sense theorem \2.2\ remains true for 'u + = oo'. 
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2.3. The jump function for a short Lin orbit segment 

According to theorem 12.21 for given a = (a~,a + ), u, A and oj = (uj~,uj + ), there is a 
unique short Lin orbit segment X = Note that X depends in particular on 

<jj, which is not reflected in our notation so far. To emphasize this dependence from we 
now use the notation X u , and similarly and v^. We define the jump function H as 

S(u>, a, u, A) := x~(a, u, A)(0) — x^(a, u, A)(0). (26) 

Using that xj + (a, u, A)(t) = q~^ + (u, A)(t) +vZ^ + (a, u, A)(t), we can write H in the form 

S(U7, a, w, A) = A) + a, u, A), (27) 

where 

A) := A)(0) - g+(u, A)(0), 

(28) 

£(u,a,u,\) := v~(a,u, A)(0) - v+(a,u, A)(0). 

Recall that (^,t>i) is the solution of the fixed point equation ( 1181) and, hence, solves 
the boundary value problem (©, (J„), (B~), (B+)). 

The term reflects the intersection of the stable and unstable manifolds of 7~ 
and 7 + respectively. We present examples for suitable choices of in section HI Here 
we focus on estimates of £. 

In order to establish those estimates, we impose some assumptions on the leading 
eigenvalues of 7~ and 7 + . Let fij denote the leading stable eigenvalue or the leading 
stable Floquet exponent of 7" depending on whether 7" is an equilibrium point or a 
periodic orbit with minimal period T~ > 0. Similarly, let /x^ be the leading unstable 
eigenvalue or the leading unstable Floquet exponent of 7 + . We assume the following: 

Hypothesis 2.1. are real and simple. 

Further, for the sake of simplicity, we also assume 
Hypothesis 2.2. dimZ = 1. 

Let Z = span {z}, \z\ = l.Then, since (eZ, 
f (w, a, u, A) = (z, a, «, A))«. 
Further we assume 

Hypothesis 2.3. The direct sum decomposition ([3]) is orthogonal with respect to the 
used scalar product (-,•)• 

Lemma 2.10. Let a, u, A, u) be in accordance with theorem \2.2\ and let hupotheses \2.1\ - 
\KEhold. Then, 

£(u>,a,u 1 \)=0(\\a\\). 

The O(-) limit holds for \\a\\ — > uniformly in u, X,uj. 
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Proof. According to the definition of £ and hypothesis 12.31 we find that 
(z,£(u>,a,u,\)) = (z,v~(a,u,X)(0) - v+(a, u, A)(0)) 
= (z, (id-P-(u,A)(0))T;-(a,«,A)(0)> 

-( 2 ,(id-P + ( M ,A)(0))^(a, M ,A)(0)). 
Since i> satisfies the fixed point equation (fT8l) . according to (fTTj) we find that 

(z, f (w, a, u, A)) = ($-(0, -u~) T z, (id - P"(u, A)(-ar))a-) 
- ($+(0, u + ) T z, (id - P+(u, X)(uj + ))a + ) 

+ (z,[ ^-(0,r) {id- p-(r))h-(r,vZ(r),u,X)dr) (29) 

+ (z, r $ + (0,r)(id-P + (r))/ i + (r,^(r),n,A)rfr). 
Jo 

First note that 

($ + (0, oj+fz, (id - P+(u, A)(^ + ))a + ) = ($+(0, c^+) T (id - P + (0)fz, a + ). 

Further, $ + (0, -) T is a solution of the adjoint of the variational equation along 7 + . 
Exponential dichotomy /trichotomy of this equation yields that, uniformly in u, A,u; + , 

($+(0, c^+) T (id - P + (0)) T z } a + ) = 0(\a+\). 

Similar arguments apply to ($~(0, — u~) T z, (id — P~(u, \)(—u~))a~). 

Standard results from Lin's method (cf. [TUl [T8l [T9] ) imply that the integral terms 
in (|29|) are also C(|a _ |) or C(|a + |) uniformly in u, A,w, respectively. Note that the 
arguments in pUlEE], where 7 are always equilibria, apply also in the present situation. 
These arguments are mainly based on the exponential dichotomy of the variational 
equation along ^ and the structure of h. □ 

Corollary 2.11. Let a + ,u, X,u> + be in accordance with corollary \2.S\ and let 
hypotheses \2. 1\W^ hold. Then, 

(z,£{u + ,a + ,u,\)) = -(<5> + (0,u + ) T z,{id- P + (u,\){u + ))a + ) +o{\a + \). 
The o(-) limit holds for \a + \ —* uniformly in u, X,u + . 

Proof. As in the proof of corollary 12. 81 we set a~ = 0. Then estimates in [191 lemma 3.20] 
provide the the corresponding estimate of the integral terms in ( J29l) . □ 

3. Joining two short Lin orbit segments 

Let 7/ U qi U 7 and 7 U q r U 7 r be consecutive short heteroclinic chain segments. The 
objective of this section is to join the related short Lin orbit segments Xi = (xf,x^~) 
and X r = (x~,x^) to a long Lin orbit segment. Here we focus on the case where 7 is 
a hyperbolic periodic orbit with minimal period T > 0. We use the same notation as 
in section [2] with an additional subscript T or 'r' referring to the left short Lin orbit 
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segments X\ or right short Lin orbit segments X r , respectively. However, for convenience 
we use the short notation uj~ = uf and uj + = uj+- 

In the construction the transition time r from Hi to S r plays a major role. In the 
present context, r is directly related to the number v of rotations of the long Lin orbit 
segment along 7. 

Theorem 3.1. Fix cu~, u + > 0. There are constants c,N>0 such that for all |A| < c, 
u = (ui, u r ), \\u\\ < c, a~, a + G W 1 , aj~ | , |a+| < c, and /or all i/ 6 N fl (iV, 00), i/iere 
a transition time r and a unique triple x = (xi,x m ,x r ), x(-) = x(v,a7,a^,u,\)(-), of 
solutions of (TJP that satisfy 

(J) xi(0),x m (0) e T,i, x m (r),x r (0) e S r , xi(0) -x m (0) e Z h x m (r) - x r (0) e Z r , 
{Bi) (id-Pf(«,,A)(-w-)) (a:,(-w-)-gf(u { ,A)(-u;-)-a-)=0, 
(S r ) (^-P+K,A)(cu+))(x r (cu+) -q+(u r ,X)(uj + ) -a+) = 0. 

Figure H] visualizes the statement of the theorem. 



im 




Figure 4. Sketch of a long Lin orbit segment (xi,x m ,x r ) near a long heteroclinic 
chain 7; U (7/ U 7 U g r U 7 r . 

We perform the proof of theorem 13.11 in two steps, see sections 13.11 and 13.21 First 
we study the flow along 7 by means of a Poincare map II : S 7 — > S 7 , where S 7 is an 
appropriate Poincare section. More precisely, we show that there are H-orbit segments y 
satisfying certain boundary conditions in E 7 . To that end, we employ a similar technique 
as used in the theory of Lin's method for discrete dynamical systems [10]. We denote 
the /-orbit that is the suspension of the Il-orbit y by x(y). Then we couple x+ and x(y) 
and simultaneously x(y) and x~ . The partial orbit x m is composed of xf, x(y) and x~ . 
The (in this context prescribed) times and u~ and the duration of x(y) add up to 
the transition time r. Further, we have x\ = xj~ and x r = x^. 

The statement of theorem 13.11 remains true for uo~ = uj + = 00 in the following 
sense: 

Corollary 3.2. There are constants c,N>0 such that for all |A| < c, u = (ui,u r ), 
\ui\ , \u r \ < c, and for all v G N fl (N,oo), there is a unique triple x = (xi,x m ,x r ), 
x(-) = x(v, u, A)(-) ; of solutions of such that for some transition time r 

(J) xi(0),x m (0) G H h x m (r),x r (0) G S r; x t (0) - x m (0) G Z h x m (r) - x r (0) G Z r , 
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(B) x l ($)eW u ( ll ), x r (0)eW s ( lr ). 

For fixed v the solution x depends smoothly on (u, A) . 

The proof of corollary 13.21 will be given at the end of section 13.21 

Definition 3.3. Ifuj~,uj + < oo we call the triple (xi,x m ,x r ) a long Lin orbit segment, 
and in case oj~ = uj + = oo we call (xi, x m , x r ) a heteroclinic Lin orbit connecting 7/ and 
7 r . If ji = j r , we call (xi,x m ,x r ) a homoclinic Lin orbit connecting 72 to itself. 

3.1. The flow near 7 

Let E 7 be a Poincare section of 7. We consider the discrete dynamical system defined 
by the Poincare map II : E T x W n — > £ 7 near p 7 := 7 (7 S 7 : 

y(n + l) = n(y(n),A). (30) 

The intersection points of qf(ui,X) and q~(u r ,X) with the Poincare section S 7 
define solutions q^(ui, X)(n), q^(u r , X)(n) of ( 1301) lying in the stable/unstable manifold 
of the hyperbolic Il-equilibrium p T . Let Ky be the (n — l)-dimensional subspace of M n 
such that 

The variational equation along has an exponential dichotomy on 7^ and we denote 
the corresponding projections by Q + (ui, X) and Q~(u r , A). Note again that the images 
of Q are well-determined: 

imQ+K A)(0) = T ?d+(u!iA)(0) V^ :A (p 7 ) C T qt(ui>xm WZ(j), 

imQ-K, A)(0) = T q - {urM0) W^ x (p,) C T q - {urtXm WX(j), 

where x denotes the (un)stable manifold of the mapping II = n(A) and we use the 
short notation at A = 0. However, there is some freedom in choosing the kernels of 

Q^, which allows us to use the ideas from Lin's method for discrete dynamical systems 
in the following, and also allows us to couple the solution of the discrete system with 
the solutions of the continuous system, cf. (149!) and (150]) below. 

Lemma 3.4. There are constants c, c > and N G N such that for all \X\ < c, 

u = (ui, u r ), U\ G Ui, u r G U r , \ui\ ,\u r \ < c, v > N and b = b~) G F 7 x Y^, \\b\\ < c, 
there is a unique solution y = y(b,ui,u r , A) of ( TffOj) that satisfies 

(B) Q+( Ul ,X)(0)(y(b,u,X)(0)-qj(u h X)(0)-b + )=0, 
Q~(u r , A)(0) (y{b, u, X)(u) - qj(u r , A)(0) - 6") = 0. 

Again we suppress the dependence of y on v from our notation. Note that this 
lemma is a discrete version of the existence and uniqueness result on Shilnikov data, cf. 
[SJ. However, also in view of its application in the following section |3~2| we consider a 
reformulation by using small perturbations of q\{n) and q^{n), similar to section [2j For 
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given i/6Nwe define v + := |_f J, the integer part of z^, and z/~ := v — u + . Further, let 
us think of y as being composed of two partial orbits as follows 

() _{y + (n), ne[o,u+)nn 
m) ~ \ y-(n-v), ne[u+,u]nN [ ' 

with the additional demand that 

y +(v+) = y-(-u~). (33) 

We write 

y + {n) = q^(ui, A)(n) + w + (n) and y~(n) = qj(u r , X)(n) + w~{n). (34) 
If y ± solve ( J30l) . then w (•) satisfy the following difference equations: 

w~{n + 1) = DiYi(q^{u r) A)(n), \)w~{n) + h~(n, w~ , u r , A), 
w +(n + 1) = D^q+im, A)(n), A)w + (n) + /i + (n, u,, A), 

where 



(35) 



h^n, w, u, A) := U(qf(u, A)(n) + iy, A) - U(qf(u, A)(n), A) 
- A)(n), A)w. 

For z/ G N let denote the space of functions {0, . . . , u} —> Y y , and let S_ u denote 
the space of functions {— v, . . . , 0} — > Y 7 . For given u + and v~ we write v := (z/ + , 
and we define the space 

W v :=S u+ x S- u -. 

Then lemma l3~4l follows from 

Lemma 3.5. There are constants e, c, c and JVeN such that for |A| < c, u :— (ui, u r ) G 
JJ\ x C/ r> wift ||u|| < c, v > N and b = (b + ,b~) £ 7 7 x Y 1 , ||fe|| < c, there is a unique 
pair w v = (w+,w~) G W v , w v = w u (b,u, X), of solutions of l[35\) in an e-neighborhood 
of G W v such that 

(B w ) Q + {u h X)(0)(w+(b,u, A)(0) - b + ) = 0, Q-(w r , A)(0)(^(6,w, A)(0) - 6") = 0, 
(C) w~(b,u, A)(-z/~) - w+(b,u, A)(z/ + ) = gj"(u { , A)(i/ + ) - g d "(u r , A)(-z/~). 

Proof. To some extent the arguments run parallel to those used in section 12.21 Here we 
only give a sketch of the proof; for more details we refer to [18] . 
First we consider the inhomogeneous equations 

w -(n + 1) = D x n{q^{u r , A)(n), A)uT(n) + 
w + (n + 1) = £>ill(g+(uf, A)(n), A)u> + (n) + 0+(n), 
with boundary conditions 

(B w ) Q + (u h A)(0)(w+(0) - b + ) = 0, Q-( Wr , A)(0)(uT(0) - 6") = 0, 

(2^) (id - Q+( U ,, A)(i/+))w+(i/+) = /?+, (id - Q"(u n \)(- v -)) w -(- v -) = /?", 
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for given f3 + G im (id — Q + {ui,\){u + )) and /3~ G im (id — Q~(u r , \){—v~)). We write 
(3 := (/3 + ,j3~). Similar to the proof of lemma |2~^ we find that the boundary value 
problem (fl3EJ),(-B&), (.£?/?)) has a unique solution w G W v , w = w(g,b, (3,u, A). 
Next we replace the boundary condition {Bp) by 

(B d ) w+(v+) - w-(-v-) =d, d G Y T 

Indeed there is a (3 = (3{d) such that w{g,b,d,u, A) := w{g,b, fl(d),u, X) is the unique 
solution of the boundary value problem (flHHl), (•£*,), The argument for this fact 
runs parallel to the corresponding construction in [19] or [TO] . 

Further, similar to the proof of lemma [2761 we consider a fixed point equation whose 
solutions also satisfy the coupling condition (C): For that we define 

d„(u, A) := gj~(u,, \){u + ) - qj(u r , X)(-v~). (37) 

Finally, we consider the fixed point equation 

w — w(li.(w, u, A), b, d v , u, A). (38) 

Here, is the discrete pendant of the Nemytskii operator defined in (TTT1) . Similar to 
the procedure in section 12721 one proves that (1381) has a unique fixed point. □ 

We define functions £y 7 xr 7 (0,c) x B UlxUr (0, c) x 5 Rm (0,c) -> Sy 7 (0,e) 

6+(6,«,A):= (id-g + (^,A)(0)) W +(6, M ,A)(0), 
61(6, u, A) := (id-g-(n r ,A)(0))^(6,n,A)(0). 

First note that 6^ depend smoothly on (6, u, A), and note further that 6jjj depend also 
on v. The 'size' of these functions is closely related to the 'size' of the jumps and £ r . 
In accordance with hypothesis 12. II we assume 

Hypothesis 3.1. The leading stable and unstable eigenvalues [i s l u of p 7 are real and 
simple. 

Note that here \i s l u denote the eigenvalues of p 1 (and not the Floquet exponents of 
7 as in section [2]). 

Hypothesis 3.2. q^{ui, A)(-) and q^{u r , A)(-) approach p 7 along the leading stable and 
unstable direction, respectively. 

Hypothesis 3.3. 6^(6, u, A) and 6j(6, u, A) are not in the strong stable subspace of the 
adjoint of the variational equation along g + (u/,A)(-) and q~(u r , A)(-), respectively, for 

71 = 0. 

Lemma 3.6. Assume hypotheses \3.1U3.3l Further, let the assumptions of lemma \3.5\ 
hold. There are functions c s ^ u = c S// "(6, u, A) such that 

\b+{b 1 u 1 \)\ = c u {b,u,\){^y u + o{\vr u ), m 

\bi{b,u,\)\ = c s {b,u,\){ f ir+o(\vr)- 

There is a constant C > such that 

\D 1 b+{b,u,X)\<C\ti u \- v , | J D 1 6j(6,w,A)| <C|/iT- (41) 
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The functions c u l s are smooth and c"/ s (0,0,0) 7^ 0. The o(-)-terms are valid for v 
tending to infinity. 

The estimates in [8] (applied to discrete systems) already provide that \b]_(b, u, A) | = 
0(\fi u \~"), cf. also corollary 13.71 but they do not give information about the leading 
term. However, this information is important for the jump estimates and consequently 
for the construction of bifurcation equations. Note that the information about leading 
terms of the derivatives in (14TI) is not needed for our purposes here, but can be computed 
in a similar manner as in [TU1 fl~9] . 

Proof. With b~\_ = u, A) := b~\_(b, u, \)/\b~\_(b, u, A)| we can write 

|6+(6,«,A)| = (St,6t(6,«,A)). (42) 

Note that (w~,w+) solves (135|) . So, applying the variation of constants formula to 
( 1361) and replacing there g^ by h finally provides 

b+(b,u,X) =^ + (0,z/+)/3+ 

v+ 

- ^ + (0, m)(id - Q + (0))h + (m - 1, w + (m - l),u h A) . 

m=l 

V V ' 



= : S 

Here \I/ + (-, •) is the transition matrix of the homogeneous equation of ( l36|) . Note that 
ty + depends on ui and A. Further, [3 + is defined by the boundary condition {Bp). 
Replacing b~]_ in the scalar product (142]) yields 

(S+ b+(b, u, A)) = (S+ M> + (0, v + ){3 + ) - (S+ S) 

= (* + (0,^) T (id-Q+(0)) T 6+,/3 + > - (S+5). 

Considerations similar to those in [TU1 [T9] show that the leading-order term of 
&+(&,«, A) will be determined by (&+ *+(0, or (*+(0, i/ + ) T (id-Q+(0)) T &+/?+>, 

respectively. Note in this respect that, due to the coupling condition (Bp), the quantity 
f3 + depends on v. 

Computations in [lOl [19] show that under hypothesis 13.31 

*+(0, „+f(id - Q + (0)) T b+ = rj + (b, u, A) (^ u y u+ + o(\^r + ), (44) 

where t] + (b,u, A) 7^ is a certain eigenvector of (DiH(p\, A)" X ) T belonging to (fi u )~ l . 
Next we consider f3 + . Combining (Bp), (Bd) and (137|) yields 

f3 + ~ F = Qd(u r , A)(-0 - q + d (ui, \)(v + ) 

-Q + (u + )w + (u + )+Q~(-u~)w~(-u~). 

We define projections Q(v) = Q(u, X)(u) by, cf. [10] for their existence, 

imQ(u) = im(id — Q + (u + )) and ker Q(v) — im(id — Q~(— 
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Applying Q(u) to (HBT) yields 

P + = Q(vM(u r ,\)(-v-)-q+(u h \)(v + )) 

- Q(v)Q + (u + )w + (u + ) + Q{v)Q-{-v-)w-{-v~). 

In [10] it has been shown that the leading-order term of the right-hand side of 
(Pf6l) is determined by the first addend, and estimates given there reveal that under 
hypothesis 13.21 we have that 

p + = V *(u,\){f)- v - + o(Kr~), (47) 

where i] u (ui, A) ^ is a certain eigenvector of DiU(p\, A) belonging to /x". 
Combining (j4"3"l) , (144|) and (J47J) provides 

(6+ &+(&, u, A)) = (t7 + (6, u, A), r/>, A))(/i u )^ + o(|^r~)- 

Here we also used that (bj_,S) = o(\/i u \~ u ); we refer again to [Uil EH] for details of the 
necessary computations. 

Finally, from the definition of rj + and r\ u it follows that 

c u (b, u, A) := ( V + (b, u, X),r] u (u, A)) ^ 0. (48) 

More precisely, linear algebra shows that r) + and r\ u cannot be orthogonal. Further, due 
to hypotheses 13.21 and 13.31 both t] + and t] u are different from zero. 
Similar computations yield the statement on 

The smoothness of c u l s follows from the smoothness of r] + and r] u . The estimates 
of the derivatives follow immediately from the considerations in [8]. □ 

The following is an immediate consequence of lemma 13.61 

Corollary 3.7. Let the assumptions of lemma \3.6\ hold. Then there is a constant C 
such that for all (b,u,X) G I?y 7 xY 7 (0,c) x B UlxUr (0,c) x £ R m(0,c) 

\b+(b,u,\)\ < c(ii u y\ \bi(b,u,\)\ < c(fi s y u . 

3.2. The coupling within E 7 

Now we have all the ingredients to couple two pairs of solutions (x~[,xf) and (x~,x+) 
of the continuous system with a solution y of the discrete system, effectively combining 
the solutions and y into one solution. 

We fix the times and uo~ at sufficiently large values by using a fixed Poincare 
section S 7 and then switch to the discrete dynamical system to describe the dynamics 
near the periodic orbit. To reflect the nature of this setting we rename cj ; + as Q + and 
u~ as Or. We choose such that q + {u h A)(f2 + ), q~(u r , A)(— Q~) G S 7 . 

In our analysis we consider xfj r as perturbations of qfj r : 

X t/r = 9{%( u I/r> A ) + v ?/r( u Vr, X )- 
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Further, we represent, as in the previous section, y as a couple of solutions (y + ,y~), cf. 
(I32p . fl33l) . which are written in the form fl34l) . The coupling is performed by searching 
for xi, x r and y such that 

x+(n + ) = y (o), x ;{-n-)=y{v), 

or equivalently in terms of the perturbations 

(C) vf(a h u h X)(Vt + ) = w + (b,u h A)(0), v~(a r ,u r , A)(-fT) = w~(b,u r , A)(0). 

The orbit x m is the suspension of the orbit y (see theorem 13. ip . Note that the transition 
time r is essentially determined by the 'length z/' of the orbit y. 

For the actual coupling analysis inside the Poincare section E 7 , we have to impose 
a technical assumption. 

Hypothesis 3.4. All solutions of (111) in a sufficiently small neighborhood of ql have 
the same transition time Q + from to £ 7 . Similarly, all solutions in a sufficiently small 
neighborhood of q~ have the same transition time Q~ from X 7 to E r . 

This can be achieved simultaneously by a scaling of the vector field in a tubular 
neighborhood along /q~ . Note that this scaling does not influence any of the previous 
results; see [H] for a similar computation. 

Hypothesis 13.41 guarantees that the points qf(ui, A)(f2 + ) + v^iuu A)(f2 + ) and 
q~(u r , A)(— Q~) + v~(u r , A)(— Q~) are both in E 7 . Further, this hypothesis allows to 
determine, cf. also (@J and ([5]): 

kerQ+KA)(0) = A)(fi + , 0)(Wf © Z,) 

= keri>+( Uj ,A)(fi + ), 

49 

kerQ-( Mr ,A)(0) = $-(u r , A) (IT, 0)(V7" © Z r ) 
= kerP-(u r ,A)(fi-), 

and (l3Tj) implies 

imQ+(«,,A)(0) cimP+(u h \)(n+), 

(50) 

im<5~(u r , A)(0) C rmP~(u r ,\)(—Q~). 

An immediate consequence of fj49l and fl50l) is (cf. also the explanations following 
lemma 12. 3j) : 

Lemma 3.8. 27ie restriction Pj~ L (ui, A)(f2 + ) |i m Q+( U! ,A)(o) ac ^ s as an isomorphism 
miQ + (u h A)(0) — > imF^Mi, A)(f2 + ). Moreover, for all v + G K, we have 

P+(u h X)(tt + )Q+(u h X)(0)v+ = P+{ Ul , A)(fi>+, 

(id-Q + (u h X)(0))v + = (id-P+( Ul ,X)(Q+))v+. 

Similarly p- r (u r ,\)(-Q-)\ imQ - {Ur:X )(o) : im<3 _ (u r , A)(0) -> imP" r (n ,A)(fi _ ) an 
isomorphism, and for all v~ aY 1 we have 

P- r (u r ,X)(-Q')Q-(u r ,X)(0)v- ^ P- r (u r ,X)(-n-)v~, 

(id-Q-{u r ,X)(0))v- = (id- p-(u r ,X)(-n-))v-. 



Heteroclinic chains involving periodic orbits 24 

Further, we use the notation 

U := Ui x U r , U 3 u = (ui,u r ). 

The following lemma is a reformulation of theorem 13. II in terms of the perturbances v l , 
vf and w ± . 

Lemma 3.9. Fix u> + , u~ > 0. There are constants Q~ , Q + > 0, N e N, c > such that 
for all |A| < c, u G U, \\u\\ < c, v > N ', and for given sufficiently small a ; ~,a+ £ W 1 
there are b 6 7 7 x Ky and af,a~ £ M n siic/i t/iat 

(Bf) (id-Pf(ui,\)(-uj~)(vf(ai,Ui,\)(-uj~)-al~) = 0, 

(zrf-P r + ( Mr ,A)(^+)«(a r , Mr ,A)(^ + )-a+) =0, 
(C) vf(a h u h A)(fi + ) = w + (b, u h A)(0), u~(a r , u n A)(-fi~) = itr(&, w r , A)(0). 

Proof. We show that (a ; + , a~ , b + , bj are uniquely determined in 

A UjA := im (id - x im (id - p-(-JT)) x imQ+(0) x imQ _ (0). 

Note that all projections appearing in the definition of A UjA depend on (ui, A) or (u r , A), 
respectively. Throughout the proof we suppose that (a ; + , a~ , b + , b) belongs to A U;A . 

We construct af, a~ and b by solving an appropriate fixed point equation. 

As a consequence of hypothesis 13.41 we know that vf(ai, ui, A)(f2 + ) £ Ky. In what 
follows we suppress the dependence on a;, uj, A and b from our notation. Therefore, in 
accordance with lemma 13.81 we find: 

= Q+(p)v+(n+) + (id - q+(o))v+(q+) 

= [unQmT 1 P^ + Wm + (id - P+(fi + h + (fi+) 

= K(^ + )|imQ + (0))" 1 < ; +a+. 

On the other hand, in accordance with lemma 13.51 

w +(o) = g+(o)w+(o) + (id - g + (o))w + (o) = b + + &+. 

Hence, vf(ai,ui, A)(f2 + ) = w + (b,ui, A)(0) if and only if 

b + = (P+(n + )\ imQ+{0) y\l l =:al l and a+ = 6+ (51) 
In a similar way we find that v~(a r ,u r , A)(— fi~) = w~(b,u r , A)(0) if and only if 

6~ = (P~ r (-fT) | im Q-(o)) 1 aj )r =: aj ir and = &J. (52) 

Altogether, for fixed itj, u r and A equations (151j) and (1521) are equivalent to the fixed 
point equation 

(a^,a7,6 + ,6~) = (&;£(&), &2(6), al/aj), a!, r K)) 

(53) 

=: £«,a ((a/ - , a", 6 + , 6"), (a,~, a+)) , 
where we consider as a mapping 

g u , : A U , A x (IR n x W 1 ) -> A„ A . 
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To solve the fixed point equation (1551) we apply the Banach fixed point theorem. 
First we show that there is a G u ,x(', a t)) -invariant closed ball B(0, e) C A u A . Then 
we prove that Gu,x('i ( a ^: a t)) * s a contraction on B(0,e). 

Let c/(iT) and c r (K) be the constants according to (1231) . The subscripts / and r 
refer to v\ and w r , respectively. With the constant c related to w, cf. lemma 13.51 we 
define 

e = e(K) :=mm{c l (K),c r (K),c}. 

Now we fix some sufficiently large K. Then the estimates given in ( 1241) and corollary 13. 71 
provide that Q + , Q~ and N can be chosen so large that for all v > N and for all 
|a^|,|a+| < e(K) the mapping G u ,x(-, {&i ■> a t)) leaves the closed ball B& uX (Q,e{K)) 
invariant. This remains true also for all larger Q + , Q~ and N. 

Due to (125|) and (jH|) the mapping G u ,x('j { a T i a t)) ls a ^ so a contraction on 
-Ba u a (0, €(K)) (with increased and N, if necessary). □ 

Corollary 3.10. Let the assumptions of lemma \3.9[ hold. Then (a^~,a~,b + ,b~) depend 
smoothly on (a7 , a+ ,u, A) . 

Proof. For fixed (u, A) the smooth dependence on (a ; ~,a+) follows by applying the 
implicit function theorem at a solution of (153]) . 

To prove the smooth dependence on (u, A) we redo the proof of lemma 13T91 to some 
extent. This time, however, we decompose and w + (0) by means of Q + (0,0)(0) 

instead of Q + (ui, A)(0). Similarly, we decompose v~(—Q~) and w~(0) by means of 
Q~(0, 0)(0). In this way we get a fixed point equation in 

A := im(id - Q + (0)) x im(id - Q~(0)) x imQ+(0) x imQ _ (0), 

where all projections are considered at (ui,X) = (0,0) or (u r ,X) = (0,0), respectively. 
Note that A does not depend on (it, A), and there is a (u, A)-dependent isomorphism 
acting between A u A and A. This leads to a fixed point equation, similar to (1531) . defined 
by a mapping 

G : A x (R n x W n ) xWxl^A. 
Exploiting this fixed point equation yields the corollary. □ 

Proof of corollary \3.2i The statement of corollary 13.21 follows immediately from 
lemma [23] with aT = a+ = 0, see also the proof of corollary 12.81 and remark [231 □ 

3.3. Jump estimates 

Let the conditions of theorem 13.11 hold, and let (xi,x m ,x r ) denote the long Lin orbit 
segment. According to ( 1261) we define: 

Ei(u,af,a+,u,X) := x t (u, of, a+, u, A)(0) - x m (y, af, a+, u, A)(0), 
H r (z/, a;~, a+, m, A) := x m {v, aT", a+, u, A)(t) — x r (z/, aT", a+, w, A)(0). 
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We now consider exemplarily the jump within X/ more closely. For that purpose we 
write Ei in the form, cf. fl27|) and (T28|) . 

E t (y, aj, a+, u, A) = A) + f,(z/, a,", a+ u, A), 

where 

Cr(u h X) :=qr(u h \)(0)-q+(u h \)(0), 
£i(v,a^,a+,u,X) := ^~(a/, u, A)(0) - ^ + (aj, u, A)(0), 

with a>i = (a~[ ,at(aj , a+ ,u, A)); cf. lemma f3T9l 

Recall that we denote the leading stable Floquet multiplier of 7 by /r\ 

Lemma 3.11. Let the constants and A be in agreement with theorem \3 . 1\ and 

let v be sufficiently large. Further, we assume hupotheses \2.1\ \2.2\ and \2.3\ hold for the 
short Lin orbit segment defined by and , and we assume that the non- orbit- flip 
condition for holds, meaning that qf is not in the strong stable manifold of '7. Then 

af, a+, u, A) = 0((^m + 0(|af I). 
T7ie 0(-)-terms are valid for v — > 00 or |a[~| — > 0, respectively. 

Proof. Lemma T2.10I yields that ^(v, a^ - , a+, w, A) = 0(|aj"|) + 0(|a+|). The coupling 
condition (C), see also fl53|) . yields that a+ = b±, and from lemma 13.61 we get 
\b+\ = c u (b, u, X)(fi u )- U + o(\fi u \-"). □ 

In what follows we assume that x\ and x r approach 7; and j r , respectively. For the 
jumps of the heteroclinic Lin orbit X\ U x m U x r connecting 7^ and 7 r (via 7) we get: 

Corollary 3.12. Let the heteroclinic Lin orbit x\ U x m U x r be in agreement with 
corollary \3.2\ and let v be sufficiently large. Further we assume hypotheses \2.1[ \2.2\ 
and \2.3\ hold for the short Lin orbit segment defined by q^ + and v~t , and we assume 
that the non- orbit- flip condition for qf holds, meaning that qt is not in the strong stable 
manifold 0/7. Then there is a smooth function cf : U x W 71 — > M such that 

(z, u, A)) = cf(u, X)(^T U + o(\fJ,T v ), 
where the o(-)-term is valid for v — > 00. 

Proof. In accordance with corollary 12.111 and ( 1291) we have 

(z, &{u, u, A)) = -(^+(0, u + f(id - P+(u, A)(0)) T 2, a+) + o(|a+|). 
Using (1531 and fj40|) (in this order) yields 
(*,G(i/,u,A)> = -($+(0, W + f(id-P+(n,A)(0))^,6i)+o(|6+|) 

= -c«(6, u, \){vl»)-"(*+(0, u+) T (id - P+( W , \){0)) T z, &1) 

+ ((// u )-") 

The notations c" and 6^ are in accordance with the proof of lemma 13.61 Note that both 
b and bj_ depend on (u, A). With that we finally get 

cf(u, A) = -c u (b, u, A) ($+(0, ^ + ) T (id - P+(u, A) (0)) T ^, S+). (55) 
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The smoothness of cf follows from the representation (I55p — recall that all ingredients 



Corollary 3.13. Let the assumptions of corollary \3.1S\ hold. Additionally let n = 3. 
Then the function cf from corollary \3.12\ satisfies c"(0, 0) ^ 0. 

Proof. Consider the explicit representation (155]) of of. 

First we make clear that c"(0, 0, 0) ^ 0: For that we recall the arguments justifying 
(j48|) and note that, in the present context (n = 3 and hence dimKy = 2), hypotheses 13.21 
and 13.31 are automatically fulfilled. 

Next we consider the second term on the right-hand side of (1551) . Due to (J2J) and 
(jl]) we have (id — P l + (0,0)(0)) T z = z. Again due to n = 3, and further due to the 
definitions of b\ and Q + , cf. ((32}, and ^2) it is clear that ($+(0, u + ) T z, 6+) ^ 0. □ 

Remark 3.14. The jump E r in S r can &e treated in a similar way. 
4. Application to EtoP cycles 

In this section we apply the theory of the existence of long Lin orbits developed in the 
previous sections to EtoP cycles. We discuss bifurcations of 1-homoclinic orbits to the 
equilibrium in the neighborhood of the EtoP cycle, and we compare these results with 
the numerical results of a concrete vector field from [Til E2], which also serves as the 
main motivating example for our studies. Here we refer to homoclinic orbits to E (near 
the cycle under consideration) making only one excursion to P as 1-homoclinic orbits. 

Let the EtoP cycle consist of a hyperbolic equilibrium E, a hyperbolic periodic 
orbit P and heteroclinic orbits qi, connecting E to P, and q r connecting P to E. Then, 
in the language of the previous sections, such a EtoP cycle can be considered as an orbit 
segment E U qi U P U q r U E. The 1-homoclinic orbits to E then can be found among 
the homoclinic Lin orbits near this orbit segment. Therefore, the bifurcation equations 
for detecting 1-homoclinic orbits to E are generated by making the jump functions 
and S r equal to zero, cf. ( 1541) : 



Note that in the present context a~[ and a+ are zero, the corresponding and £ r are 
given by corollary 13.121 

In our analysis we distinguish two types of EtoP cycles. First, we consider 
codimension-one cycles characterized by a robust heteroclinic connection q% and a 
connection q r that splits up with positive speed while moving the family parameter 
A. We prove an accumulation of 1-homoclinic orbits to E near the original EtoP cycle 
in the following sense: for each sufficiently large v e N there is a \ v for which a 1- 
homoclinic orbit exists. This homoclinic orbit performs v rotations along P before 
returning to E. The \ v accumulate at A = 0, the critical parameter value for which the 
cycle exits. 



there depend smoothly on (u, A). 



□ 



Ei(i/,u, A) = 0, S r (z/, u, A) = 0. 



(56) 
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Second, we study co dimension- two cycles, where in comparison with the previous 
one, we 'merely modify' the behavior near qf. We demand that the unstable manifold 
of E and the stable manifold of P have a quadratic tangency along The parameters 
A; and A r unfold the orbits qi and q r independently. Generically, 1-homoclinic orbits 
are still codimension-one objects — hence they are expected to appear along curves in 
the parameter space. Indeed, for each v we find those orbits on a curve k v in (A;, A r )- 
space. Each curve has a turning point which tends to the critical parameter value 
(A;, A r ) = (0, 0). Here again v counts the rotations of the 1-homoclinic orbits near P. 

Our analysis is local in nature. However, the local phenomena described above are 
part of a global scenario observed numerically in several examples. In parameter space 
1-homoclinic orbits can be continued along a curve which snakes between two curves 
(which are related to a 'quadratic tangency at q{) and accumulates on a curve segment 
for which the EtoP cycle exists, an example will be introduced in the next section. 

Finally we mention that in the recent paper [5] similar phenomena have been 
discussed. 

4-1. Unfolding of a saddle-node Hopf bifurcation with global reinjection mechanism 

In this section we consider a three-dimensional model vector field that was introduced 
in pj], the numerical results concerning the EtoP cycle presented here are from [12]. The 
vector field describes the dynamics near a saddle-node Hopf bifurcation in the presence 
of a global reinjection mechanism. This type of dynamics with reinjection can be found, 
for example, in laser systems p31 [22} [23], in dynamo theory [1] and, more generally, 
near weak resonances [211 chapter 4.3.2]. The vector-field model can be written in the 
form 

x = A\x — ujy — (ax — (3y) siny? — (x 2 + y 2 )x + d(2 cosy? + A 2 ) 2 , 

y = Aiy + ujx — (ay + (3x) sin 99 — (x 2 + y 2 )y + f(2 cosy? + A 2 ) 2 , (57) 

(p = A 2 + s(x 2 + y 2 ) + 2 cosy + c(x 2 + y 2 ) 2 , 

where Ai and A 2 are the unfolding parameters of the saddle-node Hopf bifurcation. The 
parameters u, a, (3, s, c, d and / determine the type of unfolding and we keep them 
fixed throughout at 

u = 1.0, a = -1.0, (3 = 0, s = -1.0, c = 0, d = 0.01, / = nd. 

This choice corresponds to the unfolding of type A that was studied in [11], where more 
details can be found. The variable ip is 27r-periodic and global reinjection is realized by 
trajectories that connect a neighborhood of a saddle-node Hopf point with one of its 
symmetric copies. Hence, a global reinjection corresponds to a large excursion near the 
circle S 1 = {x = y = 0}. Note that this circle is not invariant because d 7^ and / 7^ 
(where rational ratios are avoided). 

As already shown in [TT], the system has a wide variety of homoclinic orbits of 
saddle-focus equilibria involving one or more global reinjections. The most interesting 
one in the present context is the homoclinic orbit hi to the saddle-focus equilibrium 
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E = (0, 0, arccos(^2/2)) with one global excursion, which accumulates on a curve 
segment in parameter space, while the orbit itself accumulates on an EtoP cycle 
connecting E and a periodic orbit P. Figure [5] (a) shows the relevant part of the 
bifurcation diagram of (1571) where the accumulation of hi takes place. Both curves t 
are continuation curves of the co dimension-one heteroclinic orbit q\ connecting E to a 
periodic orbit P. Here the codimension is characterized by a quadratic tangency of the 
unstable manifold of E and the stable manifold of P. The curve c\ is the continuation 
curve of the heteroclinic orbit q r connecting P to E. The dimensions of the unstable 
manifold of P and of the stable manifold of E add up to the space dimension. Hence, 
q r is also a codimension-one heteroclinic orbit. 

The complete heteroclinic EtoP cycle is given by E U qi U P U q r . Panels (b)-(e) 
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Figure 5. Panel (a) shows a detail of the bifurcation diagram of ([57)) in the (Ai, A2)- 
plane. Shown arc the curve h\ of a codimension-one homoclinic orbit to E, the curve 
Ci of a codimension-one EtoP connection from E to P, and the curve to of tangencics 
of a codimension-zcro EtoP connection from P back to E. Panels (b)-(e) show the 
relevant part of a time-versus-y plot of selected homoclinic orbits on h\ that illustrate 
how they take more rotations close to P as they approach the complete EtoP cycle at 
the intersection of c\ and to. Here, T is the total integration time of the computed 
orbit segments. 
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show time-versus-;?/ plots of selected homoclinic orbits along h\ that illustrate how the 
homoclinic orbit accumulates on the EtoP cycle as the bifurcation curve h± accumulates 
on the segment of c\ where the complete EtoP cycle exists. 

Our goal here is to explain the accumulation process of h\ analytically More 
precisely, we are going to show two features of h\. First, we consider a one-parameter 
family along a curve (Ai(A), A 2 (A)) somewhere in the middle between the two curves t , 
and show that the 1-homoclinic orbit hi to E accumulates at discrete points on that 
parameter line. Second, using an unfolding of the tangencies to at the intersection point 
of to with c\, we explain the shape of hi near the turning points (near the points labeled 
(b)-(e) in figure [5] (a)) during the snaking process. 

4-2. Accumulation of homoclinic orbits near EtoP cycles 

Consider a one-parameter family of ODE (JT]), that is m = 1. We assume that for 
A = there is a heteroclinic EtoP cycle consisting of a hyperbolic equilibrium E, a 
hyperbolic periodic orbit P and heteroclinic orbits q r and qi connecting P to E and 
E to P, respectively. In accordance with the notation in section section [3], we have 
E = 7; = 7 r and P = 7. 

The aim of this section is to study homoclinic bifurcations from the given 
heteroclinic EtoP cycle under some additional genericity conditions. 

Throughout we consider the system for A G (— c, c), c sufficiently small. We assume: 

(CI) dimW s (E) = k and dimW u (P) = n - k; 

(C2) W S (E) and W U (P) intersect in an isolated connecting orbit q r ; 

(C3) The extended manifolds Ua6(-c,c) W x ( P ) x W and Ua 6 (- c ,c) W x( E ) x ( A ) 
intersect transversally in M n x R; 

(C4) W S (P) and W U (E) intersect transversally along qi\ 

(C5) The leading stable Floquet multiplier /i s (A) of P is real and simple; 

(C6) qi and q r approach E and P along the leading stable/unstable directions (non- 
orbit-flip condition). 

Remark 4.1. The following one-parameter subfamily of |57| ) satisfies the conditions 
(CI) - (C6): Let k = (Ai(A), A 2 (A)) and (Ai(0), A 2 (0)) G c x . Here 61 denotes the part 
of C\ between the intersections of c\ and the curves t ■ Further we assume that k and c\ 
intersect transversally. 

In order to apply the theory that we developed in the previous sections, we introduce 
cross-sections Ej/ r of qu r . Conditions (C1)-(C4) ensure that dimUi = dimf/ r = 0. 
Therefore, corollary 13.21 tells us that for each sufficiently small A there is a unique 
homoclinic Lin orbit x = (xi,x m ,x r ) connecting E to itself. The actual 1-homoclinic 
orbits relate to solutions of the bifurcation equation (|56|) . However, because of (Cl)- 
(C4) we have dimZ; = and dimZ r = 1. This means that the (unique) homoclinic Lin 
orbit has exactly one discontinuity, and this discontinuity is located inside E r . In other 
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words, is identically zero and the bifurcation equation for 1-homoclinic orbit reduces 
to (note that no u is involved) 

E r (u, A) = 0. 

The jump function E r is defined by ( 1261) and ( 1271) 

~r(v, A) := x m (u, A)(r) - z r (r/, A)(0) = £°(A) + A). 

Condition (C3) yields that the manifolds W U (P) and W S (E) split with non- vanishing 
velocity, in other words, -D£ r °°(0) ^ 0. Hence there is a parameter transformation such 
that 

C(A) = A. (58) 
For the remaining term £ r (z/, A) we may employ corollary 13.121 

^,A)=cf(A)o*r+o(i/' , n. 

Combining these terms yields the following lemma. 

Lemma 4.2. Assume (C1)-(C6) and 155]) . Tnen the jump function H r (i/, A) can &e 
written as 

S r (r/, A) = A + <(A) (//(A))" + o (|/i s (A)D , (59) 
where cf.(X) : R — > R smooth, and the o{-)-term is valid for v — ► oo. 
A direct consequence of this lemma is the following corollary: 



Corollary 4.3. Under the assumptions of lemma 4-2 there is a constant N E N suc/i 
£/iat /or a// ^ G N, z/ > JV, inere a A„ snc/i f/iat 5 r (i/, A^) = 0. Moreover, X u tends to 
as — > oo. 

Note that the zeros of S r correspond to 1-homoclinic orbits to E. Hence, 
corollary 14.31 says that for each sufficiently large natural v there is a 1-homoclinic orbit 
with v rotations near P. Further, we see that in parameter space these orbit accumulate 
at A = — in state space they accumulate onto the original EtoP cycle. 

For more precise assertions we need to know that c£(0) 7^ 0, which is true if n = 3, 
see in the proof of corollary 13.131 Therefore we get: 

Corollary 4.4. Let n = 3. Assume further (C1)-(C6) and 

(i) If fi s > 0, then there is a monotonically increasing/decreasing (if c s r {fS) > 0/c£(0) < 
0^ sequence (A„), A„ — > 0, such that S r (r/, A„) = 0. 

(ii) If [i s < 0, then there is an alternating sequence (A„), X u — > 0, such that 
Z r (v, \ u ) = 0. 

The order in which \ v approaches for v — * 00 is given by \i v . 

A subfamily of fl57|) as described in remark 14.11 is related to (i) of corollary 14.41 
Altogether, in respect to (I57|) corollary 14.41 explains the accumulation process of h\ 
along a line k according to remark 14.11 But it neither explains the global snaking 
behavior of hi nor the local behavior of h\ near the turning points, cf. (b) - (e) in panel 
(a) of figure El which we consider in the following section. 
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4-3. Homoclinic orbits near degenerate EtoP cycles 

In this section we consider a co dimension- two EtoP cycle. For that we modify in the 
formerly introduced EtoP cycle only the heteroclinic orbit q\. Here we assume that 
along qi the unstable manifold of E and the stable manifold of P do no longer intersect 
transversally but have a quadratic tangency. We use parameters A; and A r to unfold the 
codimension-one heteroclinic orbits qi and q r , respectively, and write A = (A r , A/) £ M. 2 . 
Here A r plays the same role as A in the previous section, and A; moves the manifolds 
W U (E) and W S (P) against each other in a direction which is orthogonal to the sum of 
the tanget spaces of these manifolds. 

We consider system (TjQ) for |A| £ (— c, c), c sufficiently small. In detail we assume 
the following: We adopt the assumptions (CI), (C2) and (C6) from the previous section 
as (CI'), (C2') and (C6'), respectively. In (C3) we only replace A by A r : 

(C3') the extended manifolds U A? . e (_ c , c )W^ r Q JP) x {A r } and l)x r e(- c ,c)W? Xr JE) x {A r } 
intersect transversally in R n x JR.. 

(C4') W S (P) and W U (E) have (along qj) a quadratic tangency, and 

the extended manifolds U Xl e(-c,c)W^ h) (E) x {A ,} and U A;e( _ CiC) W / ( s 0A;) (P) x {A,} 
intersect transversally in M. n x R. 

In contrast to (C5), here we also have an assumption on the leading unstable Floquet 
multiplier 

(C5') the leading stable and unstable Floquet multipliers /i s /"(A) of P are real and 
simple. 



Remark 4.5. Consider (El). There is a parameter transformation (A r , A;) <-» (A l5 A 2 ), 
with (Ai(0), A 2 (0)) £ ci fl to, such that in the new parameters |57| ) satisfies (CP) - 
(C6>). 

The above assumptions imply the following dimensions of the involved linear 
subspaces: dimZ r = 1, dimf/ r = 0, dim. Zi = 1, dimUi = 1. Hence, the variable 
u = ui appears in the jump function and in the bifurcation equation, respectively. 

The jump function S now consists of two parts, each representing one jump: 

E r (v,u,\) 



W, u, A) 



where E r /i(u, u, A) = ^ t (u, A) + £ r /z(^, u, A). Indeed £™ depends only on A r — more 

precisely it has (after an appropriate parameter transformation) the form, see also (158j) . 

C°°(\) = V (60) 

Further depends only on A; and u, and the quadratic tangency within can be 
modeled by 

zr(u,\i) = \i-u 2 . (6i) 

Altogether this yields 
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Lemma 4.6. Assume (C1')-(C6') and KUtyl . (E2]]. Then the bifurcation equation for 
1-homoclinic orbits can be written as 

*(uu\)=( K ) + ( (mw+^w) \ 
{ , , ) { xt-u* ) + { c ?(u,\)(^(\)r+om\)r) ) °> {b2) 

where c s l u : 1 x I 2 ^ 1, are smooth and the o(-) -terms are valid for v — > oo. 




Figure 6. Solution curves of (|63|) and (|62p in the (A r , A;)-plane for fixed v. Panel 
(a) shows the solution curve k°° of the unperturbed equation (|63p . which is a simple 
parabola in the (A;, u)-plane and a ray in the (A r , A;)-plane that is covered twice as u 
is varied. Panels (b) and (c) show possible perturbations of k°°. 



Solutions of the 'unperturbed equation' 

(w) =0 (63) 

correspond to EtoP cycles near the original one. Solutions (X r ,Xi)(u) of (1631) are 
displayed in figure El panel (a). For each A; > there are two different heteroclinic orbits 
connecting E to P. Therefore, for each nonzero A £ k°° there are two different EtoP 
cycles satisfying (CI) - (C6) (along a curve intersecting k°° transversally) . So, according 
to corollary 14.31 we expect for fixed v and A; > two different 1-homoclinic orbits 
Hi(u, Xi) and H 2 {v, X{) to E with v rotations near P. Both orbits can be continued in 
parameter space. Amazingly, they are located on the same continuation curve. Indeed, 
solutions of fl62l) for fixed v are small perturbations of the solutions of fl63|) . In parameter 
space we find those solutions on perturbations k v of k°° as displayed in figure El (b) or 
(c). Both, Hi(u,Xi) and H2(v,Xi) are on the same k v , but on different branches. For 
decreasing A;, the orbits Hi(u, A;) and H 2 (i', A;) finally 'merge' in the vertex \ u of k v . 

Generically one expects a perturbation of k°° as depicted in figure [6] (c). Below we 
show that in three-dimensional state space such a parabola like curve will indeed appear 
— as the numerical computations suggest. 

Corollary 4.7. Let n — 3. Assume further (C1')-(C6') and ffiU\) . If61\) . For each 
(sufficiently large) v £ N there is a solution curve of k v = (X r ,Xi)(u) of the bifurcation 
equation ( Tfi2j) for 1-homoclinic orbits. Further, there is a unique u u (for each v) such 
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that D\i{u v ) = and DX r {u u ) ^ 0. The points \ u : = (A r , A^f/u^) accumulate at A = 0, 
and the curvature of n u in \ u tends to infinity as v tends to infinity. 

This corollary gives an explanation of the shape of hi (in figure [5] (a)) locally around 
the turning points. It also provides information about the exponential rates with which 
the turning points accumulate to c\ fl t , namely they are given by the stable Floquet 
multiplier in the A r -direction, and by the unstable Floquet multiplier in the Ardirection. 
Similar numerical studies in [5] reveal that this mechanism also occurs in other systems. 

A possible arrangement of curves k u is displayed in figure [7J This picture verifies 
the shape of the curve hi in figure [5] near the points (b) - (e) analytically. The dashed 
line in figure UJ (the co dimension-one line of the heteroclinic orbit qi) corresponds to 
the upper curve to in figure [5] Note that the vertices of the curves k u are actually not 
located on this line. The fact that in figure [5] the points (b) - (e) are seemingly on t 
is due to the large absolute value of the unstable Floquet multiplier and the resulting 
quick convergence to t . 




A,- 



Figure 7. The parabolas u, A) = for increasing values of v in the (A r , A;)-planc. 
The vertices of the parabolas approach (0,0); the order of the displacement in the 
A r -direction is given by (pi*)" and the order of the displacement in the Aj-direction is 
given by (fi 11 )^". 



Proof. Consider the bifurcation equation ( 1621) . First we note that the derivatives of £j i r 
with respect to u and A admit the same estimate as given in corollary l3.12l for £j; we refer 
to [in] or [19] for similar assertions including proofs. So, using contraction arguments, 
we can solve S(z/, u, A) = for n u := A(w, v) for sufficiently large v. With 

£>!<((), 0)^0 (64) 

follows the existence of vertices of k u . Again using (IMl) we can write n v as A; = A/(A r ). 
From that representation one easily reads off the assertion concerning the curvature. 
It remains to verify ( 1641) . Analogously to ( 155|) . we find 

c s r (u, A) = c s (b, u, A) ($;(0, -^") T (id - P r (A)(0)) T ^, Si). (65) 

Note that in the present context {n = 3) the scalar product on the right-hand side of 
(165jl is different from zero, and the quantities within the scalar product do not depend 
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on u (dimC/ r = 0). To verify that for b J recall the definition (I39p of if[. Note that Q~ 
here does not depend on u (dimU r = 0) and dimim(id — q~(\)) = 1 (n — 3). Since 
6][ = 6]]/ we get that &[[ does depend neither on u nor on b. 

So /J^c^O, 0) is different from zero if and only if D u c s (0, 0, 0) ^ 0. Similarly to (148|) . 
we have c s (6, u, A) = (r]~(b, u, X),r] s (u, A)). Actually, here rf does not depend on (b, u). 
Further, note that rf is related to the asymptotics of \I/ _ (0, — z/~)(id — Q~(0)) T fcJ_ (see 
in the proof of lemma [376]) and none of these terms depends on u. So c s = c s (u, A), and 

D 1 c s (0,0) = (r ] ~(0),D 1 r ] s (0,0)). (66) 

Roughly speaking, r] s (u,\) is related to the asymptotics of q^(u, A)(-); see again the 
proof of lemma 13.61 

First we make clear that Di7] s (0,0) is different from zero. For this we assume 
that the traces of W S (P) in both and Sp are flat (this can always be achieved by 
appropriate transformations). Then W S (P) fl £j coincides with U, and gj~(w,0)(0) = 
4> n+ (0)(u); compare hypothesis 13.41 Actually, this mapping can be considered as a 
mapping R — > R. Since n+ (O)(-) is a diffeomorphism (R 2 — > R 2 ) also the above 
considered restriction to £7 is a diffeomorphism, and 

D lQ + (0,0) (0)^0. 

In the present context we have, see |10j . 

^(n,A) = lim (Din(0, X))- n q+(u,X)(n). 

n—>oo 

From that representation we conclude that with Diq^(0, 0)(0) 7^ also 7J'i?7' s (0, 0) ^ 0. 

Finally note that r] s (u,0) £ T^l^(p) = V7 S (P) fl S P (see above, p denotes the 
Il-equilibrium P fl Tip), which is one-dimensional. Hence _Di?7 s (0,0) points in the 
same direction as 77 s (0, 0) and, because of ( 1661) . we have DiC s (0,0) ^ (see also the 
justification of PH|) ). Therefore fl64l) holds. □ 

5. Conclusions and outlook 

We adapted Lin's method to heteroclinic chains involving periodic orbits. The main 
emphasis was on the coupling of two short Lin orbits near a periodic orbit. In this 
way, we also achieved estimates of the jump functions (Lin gaps), which are essential 
for detecting actual orbits near the primary chain among the Lin orbits. 

We employed our results to study 1-homoclinic orbits to the equilibrium near a given 
EtoP cycle. In particular we gave an analytical justification of some local phenomena 
in the course of the (global) snaking behavior of the continuation curve of 1-homoclinic 
orbits. 

A complete analytical description of the snaking behavior is still an open problem; 
a global assumption on the behavior of the stable and unstable manifolds of E and P, 
similar to that used in [2] , is necessary for such an analysis. In [5j such an assumption 
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has been used for a geometric explanation of the mentioned global snaking phenomenon. 
However, also these considerations are bound to M 3 . 

Another interesting point is the more complete description of the dynamics near 
an EtoP cycle, such as existence of iV-homo clinic (or A-hetero clinic) orbits to E or P, 
periodic orbits, or shift dynamics. 
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